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Abstract. Online parameter identification is of importance, e.g., for model predictive 
control. Since the parameters have to be identified simultaneously to the process 
of the modeled system, dynamical update laws are used for state and parameter 
estimates. Most of the existing methods for infinite dimensional systems either impose 
strong assumptions on the model or cannot handle partial observations. Therefore we 
propose and analyze an online parameter identification method that is less restrictive 
concerning the underlying model and allows for partial observations and noisy data. 
The performance of our approach is illustrated by some numerical experiments. 
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1. Introduction 

Dynamical systems like ordinary differential equations or time-dependent partial 
differential equations play an important role for modeling instationary processes in 
science and technology. Such models often contain parameters that cannot be accessed 
directly and therefore must be determined from measurements, which leads to inverse 
problems. In many applications, e.g., in model predictive control, the parameter 
identihcation has to take place during the operation of the considered system. Hence 
online methods become necessary. Examples of applications range from HVAC (heating 
ventilation airconditioning) systems via battery charge estimation to aircraft dynamics, 
see e.g. [7], [9], [H]. 

In many applications we face the additional problem of having only partial and noisy 
observations of the state. Motivated by these facts, in this paper we propose an online 
identihcation method that is also applicable in case of indirect partial observations and 
takes into account noisy data. For this purpose we employ a dynamic update law for 
both the estimated parameters and the state estimate that is strongly inspired by the 
schemes from [T] and [8]. Online parameter identihcation has been extensively studied 
in the hnite dimensional setting, e.g. [6], [10] or na. The literature becomes much 
more scarce when dealing with inhnite dimensional models as arising in the context of 
partial diherential equations. We refer to the extensive literature review in [1] and [8]. 
More recent work on this topic can e.g. be found in [2]. 

The paper is organized as follows: In section 2 we state the underlying diherential 
equation with the according assumptions and dehne the online parameter identihcation 
method. In the next section the convergence analysis of the method is discussed for 
the exact data case, the case with noisy data and also the one with smooth noisy data. 
Some examples and numerical experiments illustrate the performance of the method in 
section 4. We conclude with some remarks and an outlook in section 5. 

2. Online Parameter Identification method 

In this chapter we present the underlying diherential equation and the corresponding 
assumptions. Further we introduce an online parameter identihcation method. 

Let Q, X and Z be Hilbert spaces. We consider the abstract ordinary diherential 
equation 

gt(t,x)=0 (1) 

ufit, x) + M(t))(x) = fit, x) 

u(0, x) = u^ix) 
g(0,x) = qoix) 

where C:Q x D(C')(C Q x X) —)• X, f: [0, oo) x A —)• X and the initial value for u, 
namely uq are given. The inverse problem we are interested in is to hnd the parameter 
q from given observations of the state u over time, Guit, x) = zit, x), where G: X ^ Z 
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is the observation operator and Z the observation space. For simplicity of exposition we 
consider a linear observation operator here. Most of what follows can be carried over to 
the case of nonlinear observations. 

We will denote the exact solution by q* and u*. To define an evolution system for 
identifying q* from measurements ^ we split u* in its “observed part” Ru* = G^z G 
C V and its “unobserved” part Pu* = u* — Ru* G Af{G) C 1/ by appropriate 
projections R and P. Here H C 1/W C X and V <ZVX X with the corresponding 
embedding constants C'ypx, C'v'XX) operator Gb Z —)■ X is the 

Moore-Penrose Inverse of G. Hence the projection R for the “observed” part is the 
projection on the orthogonal complement of the nullspace of G, namely R: X ^ J\f{G)-^, 
R = G^G. The orthogonal projection P is the projection on the nullspace of G, that is 
F:X ^ W(G), P = I-R. 

Assumption 2.1. For the abstract OFF we assume that 


(i) the exact solution u* exists and stays bounded, i.e. for all times t > 0 we have 
u*{t) G B,iuo) C ViC), 

where Bp{uo) = |n + rc G H H | \\v — Ruo\\y + ||tc — Puo\\y < p|; 

(a) the operator G satisfies a Lipschitz condition with respect to the second variable, 
i.e. for all times t > 0 and for all v + w E V + V 

\\G{q* ,u*{t) +v + w) - G {q* ,u*{t))\\x < Lc {\\v\\y + HwH^) (2) 

holds; 

(Hi) the operator G can be split in a part that is dependent of q and the rest: 

G{q, u) = A(u)q + B{u); 

(iv) for all u G Bpiuo) the operator A{u): Q ^ X is linear and bounded and there exists 
Ga > 0 such that 

\\A{u* + v)\\q^x < Ca{1 + \\v\\y) \/veV (3) 

or 

||A(n* + < C'a( 1 + IMU) Vn G X (4) 

or 

||A(m* + h + < Ga( 1 + ||h||p + ||h||y) WveV,veV; (5) 

(v) there exist coercive and bounded operators M\ H —)■ X and N:V —)• X i.e. 

• there exist constants cm and Cm such that for allvEV {Mv, v)j^ > cm |kl|^ 
and ||FMn||^ < Cm Iklly; 

• there exist constants and such that for all v eV {Nv, v)j^ > 
and ||PXn||j^ < Cn ll'ylly- 

Note that by continuity of the embeddings V^X,V + V^X,@jis sufficient 
for ([3]), ([5]). Conditions 1., 2. and 3. are similar to Assumptions 1 and 2 in [8]. 

Now we want to introduce our online parameter identihcation method. Online 
identihcation means that the parameter identihcation, the data collection process and 
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the operation of the system are taking place at the same time. Accurate parameter 
values are needed for making decisions while the system is in operation. Therefore 
our online parameter identihcation method includes a dynamical update law for the 
parameter and state estimates. 


qt — A{Ru* + Pu)*{Ru — Ru*) = 0 (6) 

— /??/* 

Ut + C{fi Ru* + Pu) + + uPNPu = f (7) 

\\riU — Hu^Wy 

(g,h)(0) = (go,ho) (8) 

where uq need not coincide with uq. 

The method is strongly motivated by the methods proposed by Kiigler [8] and 
by Baumeister et. ah [1]. The main difference compared to [1] is that we also 
allow for partial observations, which often occur in applications. This is also to some 
extent possible with the method from |H], however in contrast to |S] we do not assume 
monotonicity of the operator C. 

3. Convergence Analysis 

In this chapter we consider convergence of the estimator in the exact data case as well 
as in case of noisy or smooth noisy data, respectively. To do so we take a look at the 
errors between the exact solution {q*, Ru*, Pu*) and the estimated parameter q as well 
as the error in the projected states Ru and Pu that we denote by e, r and p. The error 
components 

e = q — q*, r = Ru — Ru *, p = Pu — Pu* (9) 

satisfy the following system of differential equations, where we split up the differential 
equation for the state in the “observed” and the “unobserved” part 

et- A{u* +p)*r = (10) 

n + RC{q*, u* +p)- RC{q*,u*) + RA{u* + p)e + = 0 (11) 

Pt + PC{q*, u* Pp) — PC{q*, u*) + PA{u* + p)e + uPNPu = 0 (12) 

(e,r,p)(0) = q*,R{uQ-UQ),P{uQ-Uo)). 

Here we have used the identities Ru* + Pu = u* + p and 

C{q, Ru* + Pu)-C{q*, u*)±C{q*, Ru* + Pu) = C{q*, u*+p)-C{q*, u*) +A{u*+p)e.(13) 
as well as Assumption 12.11 

3.1. Convergence with exact data 

3.1.1. Well-definedness To obtain existence and boundedness of the solutions 

according to our method O, 0, B, we hrst multiply flTOl) and CD with e and r 
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respectively and integrate with respect to time over an interval [tififli ^2 > 0 to get, 
using Assumption 12.11 


1 
2 L 


I l|2 I II ||2 

|e||Q + IkL 


r*t2 




'*1 


= / {(et,e)Q + (rt,r)x}dr 

Jti 

(RC{q*, u* + p) — RC{q*, u*) + pRM , r'j 1 dr 


I't2 


< 


-Lc\\p\\v\\r\\x + CMf^- 


'vx 


dr. 


(14) 


iv 


We see that the equation for q was designed such that the terms containing A cancel 
out. The above estimate leads us to choose /i according to 

Assumption 3.1. For all t > 0 

Mt) > 


vx 


Therewith we obtain 


rh 




{-Lc\\p\\v\\r\\x + 2Lc\\p\\y\\r\\x} dr 


rt2 


< -Fc / Ibllyllrllxf^-r < 0 . 

Jti 

This particularly implies boundedness 


Vt>0 : ||e(t)||| + ||r(t)|||<||e(0)||| + ||r(0)|||, 


and hniteness of the integral 


VT > 0 



y\\r\\xdt < 


e(0)||^ + ||r(0)||^ 
2Lc 


< oo. 


Now it remains to hnd an appropriate bound for the error of the “unobserved” part 
of the state, which can be done quite similarly. For this purpose we multiply fll2l) with 
p and use Assumption 12.11 with ([3]) as well as flT^ to gain 

[WpWx] = iPt,p)x 

= -{PC{q*,u* + p) - PC{qfiu*),p)x + {PA{u* +p)e,p)x - {uPNPu,p)x 

< \\C{q%u*) - C{q%u* +p)\\x\\p\\x + \\A{u* + p)\\Q^x\\e\\Q\\p\\x - u{PNPu,p)x 

< ^clbllylbllx + <^^(1 + lblb)l|e||QlblU - v{PN{p + Pu*),p)x 

For the second and the last term we use Assumption 12.11 the embedding inequalities 
and Young’s inequality to get 


CAllellglblU 



+ WpWx] < 




^vxx^vvxlblblblb) 
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and 


-u{PNp,p)x < -i^cnWpW 


'ClCfi 


vx 


- iy(PNPu*,p)x < pCN\\Pu*\\y\\p\\x < V ( 


vx 


So altogether we have 


~ [llpfi] < (io + C^dklle + iCi;5;,^C™))llpllA-||p|lK 


(15) 


(16) 


+ 


C'a, 


|2 

Iq 


2 ue\\Q-^Y^Hvx 

This leads us to choose y according to 

Assumption 3.2. 


(J2 (7^ 

2 VXX ||p ^.||2 


2 c 


'N 


V 


Vt > 0 : 

to obtain 


v{t) > max < y_ , 


4(Lc + C.4(l!e(«)||g + ||p(t)||,>||p(«)||x 


Cn 


ClrCfi 


IbWII 


VX 


li[llpllt]<-pfllplfe+^^''" 




Q- 


We now define V(r(t)) := V{t) = |[|b(t)||x] and r(t) : = 


rt 


dr 

dt 


cnC% 


L y{C)d^ and hence 


vxx 


y{t). Using the former estimate we get 


< - 2 Mhci 


+ 


Cl 


vxx 


N 


-'N 


\\Pu 


. 1.2 I 

^ cnCI 


vxx 


V 


\Q 


<-V(r) + %sup||Pn*(t)||^+ 


CjY t>0 


sup \\e 


„ r<2 ,, 


2 

Q- 


Here we use the fact that for any differentiable nonnegative function p : [0, T] —)■ IRfi 
and a,b > 0 and for all t G [0,T] the following implication holds: 

p{t) < —ap{t) + h ^ pit) < - + ( 7 ( 0 ) — < max | 7 ( 0 ) 

a a \ a 


So with a = 1 and b = —^— 

cnC-^ ly 
vxx— 


2 +5Lsup,>o||Pn*(t)|U weget: 


^aPt>0 IhUj llQ I 11^ \'')\\y 

Proposition 3.3. Let Assumptions \2. 1\ with 13. il and \3.^ hold and let — q*, uq — 
Uq) E Q X (y + V). Then there exists a solution {q{t),u{t)) E Q x (V + V) for allt > 0 
and the following estimates on the parameter and state errors (cf. hold. 

(i) For all t > 0.- |le(t)||^ + ||r(f)||^ < ||e(0)||| + ||r(0)|||; 


(it) For allt >0: ||p(f)||x < max ||p(0)||x, 


cl 


(l|e(0)||| + ||r(0)||^) + §^sup4>o 11^^^ 


cjvCi 


(ill) /“ Mt)\\y\\r{t)\\xdt < < 00 
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3.1.2. State convergence In this section we will show that the estimated “observed” 
state converges towards the “observed” part of the exact solution. For improving the 
state convergence we impose an additional lower bound on p as compared to Assumption 
13.11 fnote that therewith Proposition 13.31 still remains valid). 

Assumption 3.4. There exists a constant ci > 0 such that for all t > 0 


fj.{t) > max 


2L, 


c 


Cm 


IbWII 


V ’ 


Cl r 


X 


\r{t)\\x\\r{t)\\y 


2 

VX 


Theorem 3.5 (State convergence). Under Assumptions \2.1\ with /^). 13.21 and \3.4 


we 


have that ||i?(-u(t) — u" 


X 


= ||r(t)||jj. ^0 as t ^ oo. 


Proof. We hrst take a look at the “observed” state error for t 2 > ti > 0, for which we 
get from fITTD and flT^ 


f*t2 


d 


\r{t2)rx - MhWx = I j^\\r 


\dt — 


rt2 


{rt,r)xdt 






rt2 


R(^C{q*,u)-C{qfiu +p)yr^ - (^RA{u + p)e,r'^ - (fvRM 

_ _^ ^ V 


Flip 


, r ) dt 

X 


( 1 ) 


( 2 ) 


(3) 


where we have to estimate these terms appropriately. By Assumption 12.11 the second 
term ( 2 ) can be estimated by 


\{R{C{q*,u*) - C{q*,u* +p)),r)x\ < Lc\\p\\v\\r\\x- 


Similarly for term (3) we have with Assumption fl 2 .ip with 
\{RA{u* + p)e,r)x\ < \\A{u* + p)\\q^^ 


ellg ||r||^ < + sup |b(t)||x)||e||Q||r||x 

t>o 


< ^(l|e||Q + Ibllx) 

with 


La ■■= + sup |bb)llx), 

t>o 


(17) 


which is hnite by Proposition 13.31 Using Assumptions 12 .ll and 13.41 we get for term (1) 


( fiRM- 




in / X 


< 




CMibllpx - “2Lclblblbllx- 


Iv 


Altogether we have 

ll'■(fc)llt - ik(*i)& < r {-rcii!)(i)ii7ii''(*)iu + ifmmi + ii>'(‘)iit)}* 

< ^ ^(l|e(«)|lQ + ||r(()|| 5 :)<l« < C2((2 - «i) 
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with C2 := ■^(lle(O)llQ + ||r(0)||^), which follows from Proposition 13.31 Using this 
estimate we get for any f, 7 > 0 hxed 


|r(r)||^ + (||r(f)||^-||r(r)||^)}hr 


/t-7 


< 


|r(T)ll^hr + C2 {t- T)dT = 


r 


't-7 


/t-7 


/t-7 


Hr)\\xdT + C 2 — . 


Hence we have for all t, 7 > 0 that 


k(r)|||>7lk(t)|||- 


C27 


/i-7 


From flT^ and choosing p according to Assumption 13.41 we get 


[IkWIIg+lkWIiy^S 

r^2 


hence 


h- 


cm 


vx 


dt < 


rt2 
Iti 

Cm Cl 


-Lc\\p{t)\\v\\r{t)\\x + cmP 


|r(f)|| 


2 

vx 


'dt 


V 


^£2 


\r{t)\\xdt, 


V 


, „2 11^(0)11^+ IK0)||^ ^ 

\r\\x dt < -- < 00 . 

Cm Cl 


(18) 


(19) 


( 20 ) 


We want to show that hmi_j.oo ||/’(^)|| = 0. So we suppose that hmt_^oo lk(^)ll 7^ 0. 
If this is the case then there exists a sequence (ti)iGiv with f* —)■ cx) for t —)■ 00 , and 
an e > 0 such that for alH G IV ||r(fj)|||. > e. Now we select a subsequence {ti.)j^]M 
such that for all j & IN we additionally have U. — ti._^ > Because of inequality flTS]) . 


choosing 7 = ^ we have 


2 Co 


< 


|r(r)|||dr 


'iq-7 


By summing up on both sides and using — 7 > we get for all n G IV 


n 


— <E 




\r{T)\\\dT < 


k(r)||^dr< / ||r(r)|||dT 


i=i 

which gives a contradiction to 


□ 


3.1.3. Parameter convergence The proofs in this section are to some extent similar to 
those in Section 3 of [8]. Note however, that the Lemma quantifying the relation between 
state error and parameter error can be stated in a stronger manner (cf. Lemma 13.101 
below), which enables to considerably simplify the hnal convergence proof, see Theorem 
13.111 below. In order to show that the parameter error converges to zero we start with 
some preparatory results. First we prove an estimate on the norm of the “observed” 
state error. 
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Lemma 3.6. Under Assumption \ 2. 1\ with the projected state errors r = R{u — u*) 
and p = P{u — u*) satisfy the following relation for all 0 < ta < h < tc 

rtc 


\r{tc)\\x > 


RA{u*{t) + p{T))e(ta)dT\\x 


'tb 


k(4)ll 


X - A 


'tb 


r{a)\\xda\ dr - Lc \\p{T)\\y dr - Cm p{t) dr. 


'tb 


'tb 


Proof. Integrating identity (fTTll with respect to time and using (fl^ we obtain 


r{tc) - r(4) = / rt{T)dT 


'tb 


'tb 


R{C{q*, u*) - C{q*, u* + p)) - RA{u* + p)e - jaRM- 


I v 


dr 


Taking the norm we get, using the triangle inequality and the reverse triangle inequality, 
\\r{tc)\\x + lk(4)IU > \\r{tc) - r{tb)\\x 


> 


RA{u* + p)edT\\x 


\\RC{q*,u*) - RC{q*,u* + p)\\xdT 


'tb 


'tb 


\\pRM- 


\x 


dr 


'tb 


ly 


> 


RA{u* + p)edT\\x - / LcWpWydr- / CMpdr 


'tb 


'h 


where we have used Assumption 12.11 Now we have to estimate the remaining first term 
on the right hand side. With Assumption 12.11 as well as La as in ffTTjl we get 

II [ RA{u*{t) + p{T))e{T)dT\\x 


'h 


RA{u*{t) +p(r))(e(t„) + e(r) - e(4))dr||x 


> 


> 


'tb 


'tb 


'tb 


RA{u*{t) + p{T))e{ta)dT\\x 


RA{u*{t) + p{T))e{ta)dT\\x 


WMu* + p)\\Q^x\\e{T) - e{ta)\\QdT 


'tb 


'tb 


L\ 


\r{(T)\\xd(TdT 


where we used the fact that with r > ta and (lTn|l 


|e(T) - e(ta)||Q = II / A{u* + p)*rda\\Q < La \\r{a)\\xda. 


Combining everything yields the assertion. 


□ 


Consider the right hand side in the estimate of Lemma 13.61 While by Theorem 
13.51 the negative terms containing r will tend to zero as time tends to inhnity, the hrst 
(positive) term enables us to enforce parameter convergence by means of a so-called 
persistence of excitation condition. 
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Assumption 3.7 (Persistence of Excitation). There are To,£o)70)^ > 0 such that for 
all ta > t, f & dBf{0) there exists a time instance G [ta,ta + To] such that 


I'tb+'YO 


RA{u*{t) + p{T))fdT\\x > eo . 


'h 


To control the remaining terms —Lc lb('?')|| dr and —Cm dr on the right 


hand side of the estimate in Lemma 13.61 we will combine the estimate 

r*£2 


^ "" p ||2 I ||^||2 1*2 ^ Cm 

IIQ T II llxj — n 


"^2 llrllT.^ « 


T- 


’ti iTIly 

|2 


6{t) dr, 


( 21 ) 


where 


9 = p- 


r 


'vx 


r 


V 


resulting from fllOp with some link conditions 

Assumption 3.8 (Link conditions). There exist \, k, G [1,C)o), > 0 and 

C\, Ck > 0 such that for 70 as in Assumption lS-l] the following holds. 

For all t > T\ 


Cx> 


sup 


/ l b(^)lly 

0(t) 

Mr)\\v 


A \ A-1 


A-1 

A 


nf \ 


dr 


if X = l. 


if X> I 


For all t >Ti.^ 


I't+yo 


a > 


/*(^) 


K \ K-l 


0(t) 


sup 


h(^) 


dr 


if K = l. 


if K > 1 


, rG[4,t+7o] ^('^j 

Remark 3.9. Sufficient for Assumption \3.^ is the existence of some p > 0 and a 
constant Cp such that for all t > 0 


WpWv ^ C'p||r| 




and existence of constants Cint and Cint respectively Cp and Cp such that for all t > 0 
the following interpolation estimate 

Cmtibllylbllx < Ibife < (TintlbllylblU- ( 22 ) 

and also the connecting estimate of r and p 


CpP < ||r||^“^ < Cpp 

holds. This can be seen as follows. 

Since we want to estimate the integral (^§^)we first take a look at the 
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integrand. Using the definition of 9 and the stated interpolation estimate for the state 
error as well as the connecting estimate of r and fa we get 


^ \ K—1 


< h 






and so the integral is 


rt+yo 


e{T) 


r-i+70 


dr 


< 


-dr 


4nt 


A-l 

A 


1 

Qnt 


— 1 < const. 


The second integral [ ) dr] can be estimated similarly. Again using 


It \ 6{t) 

the definition of 9 and the estimates stated in the remark yields 


||p||A ^ ^ ||^||A ^ 


'y 


< 




r 


9 Cint /^ll^llx Cifit 

Therewith the integral is hounded, using Proposition \3.3\ 


X 


r‘i+70 


l b(r)|| 

9{t) 


A \ A-l 

^ ' dr 


A-l 

A 




nt+'yo 


|r(r) 


|Ap-7^\ A-l 


lx 


dr 


A-l 

A 


< const 


provided A > 

A possible choice for n and A is to take A = and k = max{l + 2}, which 


arises from Assumption\3.f 


With these assumptions we can state the next lemma. 

Lemma 3.10. Let Assumptions 12. Jl with \3.A \3.4\ \3.1\ a.nd \3.t^ hold. 

Then, for any given 7 > 0, there are e > 0, T > 0 and Ti > 0 such that for all ti > Ti 
the following holds true: 

If the parameter error ||e(ti)||Q > 7 , then there exists a t 2 E + T] such that the 

state error ||r(t 2 )||x > £■ 


Proof. We choose Tq, Eq, 70 , t> 0 according to Assumption 13.71 fix 7 > 0 arbitrarily, 
set Ti = max{t,i} and assume that ti > Ti and ||e(ti)||Q > 7 . (Here t will be chosen 
sufficiently large below.) Setting ^ we can choose fi according to Assumption 

13.71 Now we use Lemma 13^ with ta = h, tc = h + 70 , and set ^2 = A and T = Tq + 70 
(i.e. ta = ti < tb < tb + 'Jo = t 2 = tc < ti + T = ta + T) to obtain 

/• 4+70 p(f) 

lk(A)||x = lk(A+ 7o)||x > II / RA{u*+p) c^^llx ||e(A)||Q 
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k(4)ll 


X 


rtb+io rr 
Jh Jti 

> £o||e(ti)||Q - \\r{tb)\\x - L 

rh+'ro 
'tb 


r{a)\\xdadT - Lq 


rtb+io 


r^b+io 


ttb 


\\p{T)\\ydT - Cm / fi{T)dT 


'tb 


r{a)\\xdadT 


-L 


c 


\\p{T)\\ydT - Cm 


'tb 
rtb+'yo 


p,(T)dT. 


'tb 


The last three terms remain to be estimated. 


rtb+io /•'T 


\r{a)\\xdadT < 


rtb+'yo r-r 


'tb 


'ti 


Fb 


sup \\r{a)\\xdadT < 70 Tsup ||r(cr)||x- 

Iti (T>ti cr>ti 


Estimating by Holder’s inequality and using the link conditions results 


m 


ptb+yo 


||p(r)|| dr = 


rtb+yo 


'tb 


'tb 


lb(^)lly 

0 (r)^ 


9^{T)dT < C\ 


rtb+yo 


9{T)d 




and analogously for the last term 


^■4+70 


fi{T)dT = 




'tb 


'tb 


9{t )« 


-9-{T)dT < C^ 


'tb 


rtb+yo 


'tb 


9{T)dT 


Now using fl^ we can estimate the term 9{t) dr as follows 


^4+70 


9{t) dr < -;;^[||e||Q + ||r||^]*^+^° < ;;^(||e(4)||Q - ||e(4+ 7 o)||q) + ;;^lk(4)||x- 

I T\ ^ /Jl/f I H /T 


/tb Cm Cm 

At this point we utilize flTOl) . Assumptions 12.11 and [T7I as well as Proposition 13.31 

|2 


r* 6 + 7 o 


|e(4)||Q - ||e(4 + 7 o)||Q = [||e(t)||yt^+^o =-2 / (ei,e)Qdr 


r4+7o 


< 2 




l|et||Qlle||Qdr = 2 


^•*6+70 


Fb 


rtb+yo 


tb 

\\A{u* + pyr\\Q\\e\\QdT 

rtb+yo 


<2La I \\r{T)\\x\\e{T)\\QdT <2La I IkWIIxdr J||e(0)||^ + ||r(0)||^. 

Jtb Jtb '' 

Hence altogether we end up with 
\\r{t 2 )\\x > eo7 - \\r{tb)\\x - L\ioT sup ||r((T)||x 

/ 2 La 


(T>tl 


—LcC\ 


-CmC, 


V . 

2La 


h 


k(^)IU(irJ||e(0)||| + ||r(0)||^ + —||r(4 
V Cm 


rh +10 


Cm Jtb 


r{r)\\xdTJ\\e{0)\\l + ||r(0)||^ + —||r(4 


By Theorem 13.51 for t sufficiently large, 4, H, r > Ti > t the sum of all negative terms 
will be contained in the interval [—^, 0 ], so that we get 


I Ml eo7 

\r{t 2 )\\x >co-i- — = —. 


With e = ^, this implies the assertion. 


□ 
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Theorem 3.11 (Parameter convergence). Under Assumptions \2.1\ with \3AA \3.4 
3. 1\ and \3.8\ we have that 


||g(t) — q*\\Q —!• 0 as t —)■ oo . 

Proof. Contraposition in Lemma [3.101 gives the following assertion (as we have imposed 
Assumptions 12.1113.2113.4113.71 and 13.81 to hold): 

For any given 7 > 0, there are e, T, Ti > 0 such that for all p > Ti the 
following holds true: 

If for all t 2 G + T] the state error ||r(t 2 )||x < then the parameter error 
l|e(ti)||Q < 7- 

Thus, given arbitrary 7 > 0, we choose e and Ti according to Lemma 13.101 Then, by 
Theorem 13.51 there exists T 2 > Ti such that for all f > T 2 we have ||r(t)||x < Hence, 
for all ti > T 2 , the above statement yields ||e(fi )IIq<7- □ 


3.2. Convergence with noisy data 

In case noisy data are given instead of 2 : and the range of G is non-closed, the quantity 
might not be well-dehned, and even if it is well-dehned it will not depend on z^ in 
a stable manner. Thus we dehne a regularized version of the “observed” part of u* 

ui = GaZ^ 

with Ga a regularized version of G^ with regularization parameter a, dehned, e.g., by 
the Tikhonov-Philips method 

G^ = {G*G + aI)-^G* : Z Af{G)^ C X 

with G* : Z ^ X the Hilbert space adjoint oi G : X ^ Z, and a > 0 appropriately 
chosen. Additionally one might add a stabilizing term dehned by another parameter 
a = a{t) > 0, see e.g. |6]. Note that also the case cr = 0 is included in our analysis. 
As a matter of fact, it turns out that this term is not really needed. For the sake 
of completeness to some extent we will also consider the case of strictly positive a. 
The case of partially vanishing, partially positive a is not included here (but could be 
approximated by some positive a which partially gets arbitrarily small). 

Therewith, we redehne the estimators g, u by 

fi-A{ul^ +PuY{Ru-u^^) =-aq (23) 

ut + G{q, ui + Pit) + URM ^ + oPNPu = f (24) 

(g,'a)(0) = (go,Mo) (25) 

where a = a{t), p = p,{t) and u = fit) are chosen properly dependent on the noise level 
5{t) in 


m > ii3'(() - 3(()iu. 


(26) 
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3.2.1. Well-definedness For showing well-definedness we take again a look at the error 
components e = q — q*, r = Ru — Ru*, p = Pu — Pu* and the errors including the 


regularized version of the “observed” part 

= Ru — = r — df^ and d^^ = — Ru*. (27) 

Therewith the equalities 

-|- Pu = «„ + Pu u* — Ru* — Pu* = u* + d^^-\- p (28) 

hold. Then the differential equations for the errors are 

et - A{u* + di+ p)*rl = -aq (29) 

rt + RC{q*, u* + di+p)- RC{q*, u*) + RA{u* + 4 + p)e + piRM-^ = 0 (30) 

\va\\v 

Pt + PC{q*, u* + d^^-\-p) - PC{q*, u*) PA{u* + + p)e vPNPu = 0 (31) 

(e,r,p)(0) = {%-q*,R{uQ-Uo),P{uo-UQ)). (32) 


In case of noisy data we get a wellposedness result too. As in the exact data case some 
assumptions concerning the parameters p and v are required. 

Assumption 3.12. For all t > 0 


(^) 


pit) > max 

, 4Ca 


4L 


c 


Cm 


(Il4(t)||p + |b(t)||p)||r(t)|| 




Cm 


(l + ll<Wllv + lbb)lly)l|e(f)||Q||<(t)|| 


X, 


2cr(t) 


Cm 


\rim 


V 


2 ■ 

yx 


ftzj 


u(t) > max |z/, + ^ ||^. 


2CaC. 


Cn 

vVx^vxx 


illy) 


Cn 


IbWIlyj 


A lbb)ll 


X 


lbb)l 


yx 


A condition on the error between the regularized version of the “observed” part 
and the exact state is also needed, namely for all considered time instances t 


Il4(t)||y< 


Cm 


2 

yx 


(33) 


‘2Cm Ibl'-dlx 

should hold. This condition on smallness can be further accessed using the fact that 
da = GaZ^ — Gb and fl26l) . based on results of regularization theory and an appropriate 
choice of a{t) in dependence of 6{t) and z^{t), see, e.g. [3]. We now prove that q and u 
according to (|23il and (|24ll are well defined at least up to a certain time. 
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Proposition 3.13. Let Assumptions \2.1] with ^ and \d.l‘A hold and let — q*,UQ — 
Uq) E Q X iy + V). Then there exists a solution {q{t),u{t)) E Q x iy + V) of 
for all times 0 < t < T* where 

T* = minjt > 0 : ||4||t/> • (34) 

I Irllx I 

(i.e. the first time, when condition is violated) and satisfies the following error 
bounds (cf. ^)- 


(i) Case: a = O; For all 0 < t < T*: ||e(t)||Q + ||T’(t)||x ^ I|c(0)||q + ||?"(0)||x < oc; 
Case: a > O; For all 0 < t < T*: 

MWq + Mx < max {||g*||^, |le(0)||| + ||r(0)||^} < oo; 

(a) For all 0 < t < T*: 


lb(i)ll^<2 


clci~ 

^ VXX 

-- 


SUPi>o||^W*(OII| + 


CaC^~ 

VXX 

cnE 


[l|e(0)|b + ||r(0)||^],|b(0)|| 


(iii) IfT* = oo (cf anda = 0 then \\p{r)\\y\\r{T)\\xdT < -— 

Proof. 1. For proving the proposition, like in the exact data case we take a look at the 
norms of the squared errors. 


^^[Ibllg + Ibllx] = (et,e)Q + {rt,r)x 
= {A{u* + dl+pyrl,e) - {RA{u* + di+p)e,r) 


X 


( 1 ) 


(ab e)Q + {RC{q\ u*) - RCjq*, u* + di+ p),r)^^ - /i (^RM- 


( 2 ) 


(3) 


\'aWV /X 

"v* 

(4) 


Our goal is now to estimate all these terms appropriately, 
ad (1): Using the identity r = and Assumption 12. II we get 

{A{u* + di+pyri,e)Q- {RA{u* + di + p)e,ri + di)^ = - {RA{u* + di + p)e, di)^ 

= - {RA{u* + +p)e, d^)^ < C'a(1 + ||d^|b + lbllp)l|e||Q||d^||x 

ad (2): With some computations we get 


-(^b e)Q = -(7{q ± q*, e)Q = a{q*, e)^ - a(e, e)Q 
< -^l|e||J + cT\\q*\\Q\\e\\Q < -a||e||J + ^{\\q% + 


a 


Ik 


* 112 
Q 


a 

2 


e 


||2 

Wq- 


ad (3): The Lipschitz condition on C yields 

(RC{q%u‘) - RC(q\n' + di+p),r)„ < \\C{q\u') - C(g\u’ + 4 +P)lkl|l-|U 

^ rc(||4lli> + IIpIIk)II’'II-v- 
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ad (4): Using coercivity and boundedness of M fAssumption 12.1() and t < T* with T* 
as in fl34|) results in 


H RM-. 


a 


r = 


'' {RMr,r)^ + -;Aw7i^^Cr) 


< —p,CM 


'ally /X ll'ally 

1^11 ? _ 11 ^5 I 

I ''vx 


r 


oily 


X 


r 


5 II - 
ally 


+ itCm 


It/Fllx 


y^S I I _ 

ally 


< —TTCm 


'yx 


r 


<5II - 
ally 


(36) 


Inserting in the original inequality gives 

4i|||e||J+||r||y <C^(l + ||4||t.+ ||p||,>)||e 


5 II , ^ 

^+2 


* 112 
Q 


a . 
— e 


Q 




+Lci\K\\y+\\p\\y)\\r\\x-^CM 


'yx 




<5II , 
oily 


(36) 


Using Assumption 13.121 on /x we get 


^ ^ ril „l |2 I ll.„l|2 1 ^ II „*||2 


dt 2 


IQ 


+ l|r||i]< 


O' , 


Q 2 ^Q 


/iCM Ir Ilyx ^ ^ 


\r 


a I — 

5 lly 


Q 


-(l|e||^ + ||r||^)). 


Now we distinguish between the two cases cr = 0 and a > 0. For the hrst case cr = 0 we 
have 

4i[||e||| + ||r||y < 0 ^ i[||e||^ + ||r||l,] < i||le(0)||J + ||r(0)||i.l. 

For the second case cr > 0 we dehne r(t) := Jq o(f)df, V(t) := ^[||e(t)||Q + lk(^)||x] 
V(T(t)) := V(t). Differentiating V with respect to r leads to 


< 




,* 112 


i - v(t(«)). 


^V(r(Q4ll|e||J + l|r&l^ 

So we have for all t > 0 

^[l|e||J + Ikllx] <max|^||g*||^,^[||e(0)||J + ||r(0)||^]| < cx). 

2. We now consider the error for the “unobserved” part of the state. Similarly to ffT6|l 
we get 

^^[Ibllx] = {Pt,p)x 

- {PC{q*,u* + di+p)- PC{q*,u*),p)j, - {PA{u* + di+p)e,p)^-iy {PNPu,p)^ 

< LciWdiWy + \\p\\y)\\p\\x + Ca{ 1 + \\di\\y + \\p\\y)\\e\\Q\\p\\x 
UCn 


■INI 


yx 


+ y- 


~^c 


■||P«*|| 


■N 


V 


< ^(IIcIIq + C'0^x<^yyxlbllylblN) + (^c + C'A||e||Q)(||4lly + lblly)lbllx 


2 

ycAT 


Ibl 


(^2 

2 + ^ yXX ||p^*||2 


y 


2 2cAr 

yc + c^NbllQ)(lb^lb + lblly) + 


5 II , ll„ll 6 , ^ACyxxCyyx ||^|| 


y 


INI 
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PCIC‘1^ 

2 + ^ y^^ wpu 


* 112 

y* 


2 2 " 2cn 

Here we have used Assumption 12.11 and ([15]). Now we make use of Assumption 13.12l on 
V to get 

Ca. 


~i\\pm < -^ibi 


'N 


* 112 

lip 


|2 

l<3' 


We again define functions r(f) := •= |lb(f)llx • = 

^VXX 

V(t). Differentiating V with respect to r leads to 




PCn 


Ibl 


2 .112 ^^||ell2 

^ 2c^ ll^« lly+ 2 


vx 


2Cl~~ 

vxx 
CNl^it) 


/~i2 fi2 f-ii 

■^ibiipv+ "" 


* 112 


CaCI 


vxx 


-N 


^ ^ CNv{t) 


IQ- 


Using the embedding VX ^ X and the estimate for u in Assumption 13.12l gives 


^V(r) < -V(r(f)) + 




vxx 


IIPu 




* 112 
ly 


CaCIt 


+ 


vxx 


CnP 


\Q- 


From this we get the assertion. 

3. This is a consequence of inequality 

^ llp(T)lloll>'(T)|Lv(ir < /j(r) 

<5llWO)IIJ+l|r(0)||5,]. 


and Assumption 13.121 


|r(r) 


'vx 


Idfr) 


dr < 


IV 


^^[lbb)||Q + ||r(r)||^]dr 


□ 


3.2.2. State convergence As in the exact data case we introduce an additional lower 
bound on /i for proving convergence of the “observed” part of the state estimate. 


Assumption 3.14. There exists a constant ci > 0 such that for all t > 0 


fi{t) > max 


4L 


c 


Cm 
^ 4Ca 

Cm 


(Il4b)lb + lbb)lb)lbb)ll 


X 


(i + ll4b)lb + lbb)lb)lb(t)||Q||4lU,ci||r 


V 


2 ' 
VX 


Theorem 3.15 (State convergence). Under Assumptions \ 2. 1\ with ^ and \3.1J\ and if 
T* = oo (cf. and a = 0 we have that ||i7('u(f) — 'U*(f))||x = |b(^)llx —t 0 as f —)■ oo. 


Proof. The proof is quite similar to the one in the exact data case. We start with 
considering the “observed” state error for ^2 > U > 0. 
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r-t2 


\r{t 2 )rx - MhWx = I irt,r)x 

rt2 


'ti 


'^1 


< 


{RC{q\ u*) - RC{q*, u* + di+p)- RA{u* + 4 + p)e - -r-^RMri, r) dr 


I't2 


Lci\\di\\v + lblly)lkllx + C'a( 1 + IblU + lballx)||e||Q||r||x 


fiCM lb 


Wx 


11 — 
'ally 


dr 


Here we have used the identities and Assumptions 12.ll and T* = oo. Furthermore 
we will denote 

La ■= C'a(1 + sup{||4||x + IblU})- 


t>o 


Assumption 13.141 and Propostition 13.131 give us 

lbb2)iil-lbbi)ll^ 


f*t2 


< / {-^cdb^lb + lblly)lblU + ^A||e||Q||r||x “ 2Ac(|b„|b + lblly)lbllx}dT 




< 


< 








'ti 


4q + Ibllx)^^^ 


|e(0)||^ + ||r(0)||^)dr = C2b2-H) 


where we have dehned C 2 := {lb( 0 )||Q + |b( 0 )llx}- 

As in the exact data case (cf (ITSjl ) we get for any hxed t, 7 > 0 


't-7 


\r{T)fxdT > 7|bb)llx - 


For cr = 0 the proof from now on is exactly the same as in the exact data case. 


□ 


Remark 3.16. If in ^ T* < 00 we cannot expect convergence of the state error to 
zero if S > 0. However in this case the definition ofT* implies 


Il4(miy > 


CM MT*)\\h 
2Cm \\riT*)\\x 


and therefore that r{T*) is small, namely in case the interpolation inequality holds 
we even have that at time T* the “observed” state error is (up to a constant factor 
■A£j^) dg small as the error in the “observed” state, both of them in the V-norm. 

CintCM' ’ ■' 


3.2.3. Parameter convergence For proving that the estimated parameter converges to 
the exact one we again need two Lemmas. 












An online parameter identification method for time dependent PDFs 


19 


Lemma 3.17. Under Assumption \2. 1\ with ^ the projected state errors r and p satisfy 
the following relation for all 0 < ta < tb < tc: 


k(4)IU > II / RA{u*+ di+p)e{ta)dT\\x-\\r{tb)\\x - L\ 


riWxdsdr 


'tb 


• tb ta 


-La 


a\\q\\QdsdT-Lc {\\di\\v-L \\p\\v)dT - Cm hdr. 




'4 


'h 


Proof. The proof is basically the same as in the exact data case with ||d^||p + |b||y 
instead of ||p||y and || instead of ||r||x in the term with La and the additional term 
with a. □ 


The persistence of excitation condition is nearly the same as in the exact data case, 
except that we have instead of Ru*, i.e., here we have the regularized version of the 
“observed” part of the state. 

Assumption 3.18 (Persistence of Excitation). There are Tq, eo,'yo,t > 0 such that for 
all ta>t,f,E dB^(0) there exists a time instance tb G [tafia + ^o] such that 


rtb+io 


RA{u* + + pfififdrWx > Sq. 


'tb 


Also the link conditions are quite similar. With a slightly different definition of 


theta 


6 := /X- 


'vx 


\r‘ 


(37) 


oily 


and involving the error between the exact “observed” part and its regularized version 
we use the following link conditions. 

Assumption 3.19 (Link conditions). There exist A, k G [l,cx)), T^, > 0 and C^, 

Cfi > 0 such that for jq > 0 as in Assumption \S.lSi the following holds. 

For all t > T^ 


C-x>( 


sup 

^ re[t,t+'yo] 


( (Il4(r)llv- + llp(r)|lv0 

( fl(T) 

ll<(T)llv + l|p(r)||,> 


A \ A-l 


A-1 

A 


efi) 


dr 


ifX = l. 


if X> I 


For all t >Tji 


r-t+70 


a > 


efi) 


dr 


if K > 1 


hdj)_ 

, re[t,t+ 7 o] d(r) 


sup 


if K = l. 
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Furtheron we just consider the case (7 = 0. In the other case cr > 0 we cannot prove 
parameter convergence. The second lemma that is needed for parameter convergence is 
exactly the same as in the exact data case, (cf Lemma I3.10p 


Lemma 3.20. Let Assum,ptions \2. 1\ with /^). \3.1^\3.1f \ \37TR a,nd \S.19\ hold and (7 = 0. 
Then, for any given 7 > 0, there are £ > 0, T > 0 and Ti > 0 such that for all ti > Ti 
the following holds true: 

If the parameter error ||e(ti)||Q > 7 , then there exists a t 2 E + T] such that the 

state error ||r(t 2 )||x: > £• 


Proof. In case (7 = 0 Lemma 13.171 with Assumption 13.181 gives the same estimate as 
in the exact data case with the only difference that we have to replace with La 
and ||p(T)||y by (IIp(t)II^ + ||<i‘|lv). Thus with the adaptations we have made in the 
dehnition of 9 and in the link conditions 13.191 the proof obviously goes through like the 
one of Lemma 13.101 □ 


For T* = 00 we can prove parameter convergence analogously to Theorem 13.111 

Theorem 3.21 (Parameter convergence). Under Assumptions \2. 1\ with \3.1^ 3.14 
\3.18i. and \3.1tA and if T* = 00 (cf. ^), a = 0 we have that 


||g(t) — —!• 0 as f —)■ 00 . 

Proof. See exact data case. □ 

Remark 3.22. In case T* < 00 we cannot prove parameter convergence, because in the 
persistence of excitation assumption we need to have f —)■ cx). 

Since in case T* < 00 we cannot completely prove convergence for the noisy data 
case we also take a look at the smoothed noisy data case. 


3.3. Convergence with smoothed noisy data 

With smoothed noisy data we denote that is smooth with respect to time (which can 
be achieved by averaging over sufficiently large time intervals), i.e. 

\\{z^-z)M\z=\\zf{t)-zfit)\\z<m- 

That means for the error of the “observed” part of the state 

n = [ri + (ia]t = + [GaZ^ — z]t = + d^ 

with (i^ := Gaizf — Zt) + (Gq, — G^)zt + afi^Gcfjz^. Therewith the online parameter 
identification method as in Section is given by 


qt - A{u* + di+ p)*{Ru -ui) = 0 

(38) 

, X s Ru — 


Ut-\-C[q,u -\-d^-\-p')-\-frRM ^ -\-vPJ31Pu — f 

\\Ru — u^Wy 

(39) 

(g,M)(0) = {qo,Uo) 

(40) 
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Hence we can alternatively to consider the equations 

ct - A{u* + + pfirl^ = 0 

rit + di + RC{q*,u* + di+p)- RC{q\u*) + RA{u* + < + p)e + = 0(41) 

Ir oil V 

Pt + PC{q*, u* + d^^+ p) — PC{q*, u*) + PA{u* + + p)e + vPNPu = 0 

(e,r^,p)(0) = {% - qfiRu{R) - G„/(0),P(mo - Mq)) , 

that upon the replacements r i—)■ r^, w* i—)■ u* + and up to the perturbation in fj4Tl) 

are the same as (HOD - (ED. 

Here and below we set cr = 0. Since now we only deal with (and not with r, 
simultaneously as in the previous section) proofs become much more analogous to the 
exact data case. 


3.3.1. Well-definedness For proving the well-dehnedness, again some lower bounds on 
pi and p are required. 


Assumption 3.23. For all t > 0 

(^) 


pi{t) > —(Lc(lld^(t)llp + Ib(t)llp) + lldl(t)IU) 
Cm 


xK 


V. 


2 

vx 


(ll) 


p{t) > max |z/, ^ 


4(Lc + CA\\e{t)\\Q) 


2CaC. 


Cn 

vvx^vxx 


mt)\\v + \\d: 


^aWV) 


Cn 


IbWIlvj 


A lbb)ll 


X 


lbb)l 


vx 


Therewith we get a similar well-posedness result as in the previous section. The 
critical difference to Section 3.2 is that we get existence for all times. 

Proposition 3.24. Let Assumptions 12.II with ^ and \3. 2,1l hold, and let (go — q*, ho — 
Mo) G Q X (H + V). Then there exists a solution {q{t),u{t)) G Q x (H + H) of (E3), 
[4^ ) for all times t > 0 and satisfies the following error bounds (cf. ^). 

(i) For allt>t) : ||e||^ + ||ri||x < l|e(0)||| + |ba(0)||x 


(ll) Forallt>t) : |b(t)||x < max <{ |b(0)||^, sup^^o ll^w*b)ll 




CaC: 


vxx 

"2- 


(ill) bTlbWIblbc 




,5 


(l|e(0)||^ + |bi(0)|||: 


r^(0)||^ 


2Lc 


< OO. 


3.3.2. State convergence To obtain state convergence we have to replace the parameter 
choice from the exact data case in a straightforward manner with replacements |blly ^ 

WpWv + ll'^ally + Ibally + r ^ r)^. In case of smoothed noisy data we do not 

need conditions on ||d^||y and therefore we can prove state convergence to 0 as f —)■ cx). 
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Assumption 3.25. There exists a constant Ci > 0 such that for all t > 0 

fi{t) > max | —(Lc(||4|ly + \\p\\v) + \\di\\x ), Ci||r^(t)||x| 1^2° 

L^M J \Ta{^)\\y 

Theorem 3.26 (State convergence). Under Assumptions \2.1\ with \3. 2 A and \3.2^ 

we have that ||-R'u(t) — GaZ^{t)\\^ = ||r^||x —t 0 as t ^ 00 . 


3.3.3. Parameter convergence Due to the inhomogeneity in the right hand side 
of the equation for the “observed” state error (HTD the crucial estimate for parameter 
convergence becomes 


Lemma 3.27. Under Assumption \2. 1\ with the projected state errors rf^ = Ru—GaZ^ 
and p satisfy the following relation for all 0 < ta < h < tc ■ 


K{tc)\\x> 

ptc pr 


RA{u*{T)+di{T)+p{T))e{ta)dT\\x-\\ri{tb)\\x- / ||di(r)||xdr 


'4 


'4 




’ / lK(s)IIX dsdr - Lc i\\d^{r)\\v + dr - Gm /i(r) dr 

To obtain parameter convergence we use the persistence of excitation and link 
conditions, Assumptions 13.18) and ITT^ with the only slight modihcation as compared to 


r.S 112 


9 = 9 ^ = Pi- 


<vx 


\r 


aUV 


Therewith we obtain: 


Lemma 3.28. Let Assumptions \ 2. 1[ Id. id! fXTPl \3.23l and \3.2^ hold. 

Then, for any given 7 > ^ sup^^g ll^i(^)||v; there are £ > 0, T > 0 and Ti > 0 such 
that for all ti > Ti the following holds true: 

If the parameter error |le(fi)||Q > 7 , then there exists a t 2 E [ti,ti + T] such that the 
state error ||r^(t 2 )||v: > 

This allows us to conclude: 

Theorem 3.29 (Parameter convergence). Under Assumptions \2.1[ \3.18[ Id. idl Id.^dt 
and \3. d,5l we have that 

limsup \\q{t) - q*\\Q < — sup \\di{t)\\x ■ 

t^oa £0 i>0 

Proof. By contraposition in Lemma r3.28( for all 7 > ^ sup^^g ||d^(t) ||x there exists 
e > 0, T > 0, Ti > 0 such that for all ti > Ti and for all ^2 ^ [hfii + T]- ||r^(^2)||x < ^ 
implies ||e(f)||Q << 7 . 

So for given 7 > ^sup^^g ||d^(t)||x: we choose £ and Ti > 0 according to Lemma 
13.281 Then from Theorem 13.261 it follows that there exists T 2 > Ti such that for all 
f > T 2 we have ||r;^(t 2 )||x < hence by the above ||e(fi)||Q < 7 for all R > Ti. Since 7 
can be chosen arbitrarily close to ^ sup^^g ||d^(t)||x the assertion follows. 


□ 










































An online parameter identification method for time dependent PDFs 

4. Numerical experiments 


23 


4 . 1 . Identification of a coefficient in a degenerate diffusion equation 

Consider the problem of identifying q = q{x) on a domain C in the (possibly 
degenerate) parabolic initial bonndary valne problem 

ufit, x) — V {D{x)'\/u{t, x)) + q{x)u{t, x) = fit, x) in (42) 

u(t, x) = g(t, x) on dQ 
u(0, x) = uo(x) 

from measnrements of the state n on a snbdomain u C Q 

z(t, x) = Gu(t, x) = u(t, x)jaj (43) 

Here f(t) E -h^(ff), g(t) E D E L°°(H), are assnmed to be known and chosen 

snch that for q = q* a solntion u{t) = u*{t) E iif^(f2) to fl42]) exists for all times t > 0. 
With the spaces 

Q = HfiH ), X = L^n) , Z = L‘^{u) , 

where s > ^ so that Q is continnonsly embedded in L°°(ff), and the operators dehned 
by 

C{q,u) = B{u) + A{u)q, {B{u),v)x = 

this hts into the framework of the previons sections with an appropriate choice of the 
spaces V,V,VX,VX and the operators M,N, see below. Note that this formnlation 
corresponds to the standard semigronp formnlation for parabolic problems in case D > 0 
(see, e.g., m- However we do not assnme D to be positive, not even nonnegative, 
hence the monotonicity assnmption from [5] fails even if we nse the setting there with 
the problem adapted spaces V = {v E L^(f2) | ^yiDlVv E L^(f2)} with the norm 

ll'^llv' = ^11 \/\D\Xv\\‘^ 2 (^q) + j ^, H = Lfi{Vt) (with the notation V and H from 

0)- The case B < 0, often denoted as antidiffnsion, for example occnrs in certain 
models of pattern formation, see e.g. [5]. 

We hrst of all dehne the spaces V, V snch that the Lipschitz condition on C from 
Assnmption 12.11 holds: 

V = {v E Lf{Vt) I snpp(il)Vn) C ix , snpp v C u and V(T)Vn|tj) G L‘^{oj)} 

V = {v E ififil) I snpp(DVn) C \ a;, snpp n C \ ca and V(DVn|n\tj) E ififil \ a;)} 
with norms 

ll'^llv ~ II+ ||n||i 2 (Q), ||n||y = ||V(DVn|Q\tj)||L 2 (f 2 \^) + ||n||i 2 (f 7 ). 

and their smooth connterparts 

= {4> E Cf°{Q) I snpp(0) is a compact snbset of uj} 

= {fj E (7^(12) I snpp(';/^) is a compact snbset of \ w} 


{DVuY'Vvdx, A{u)q = qu , Gu = u\^ 
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which are dense in V and V, and whose sum id + is dense in Therewith the 

Lipschitz condition on C{q*, •) is obtained as follows. For any (p + ^ V + V we have 

{C{q*, u*{t) + v + w) - C{q*, u*{t)), + 'fi)x 

= / ({DViv + + ip) + q*{v + w){(j) + fi)) dx 


= / q*{v + dx + {D'Vv\uj)'^'V(j)\ujdx + / {D'Vw\n\,j)'^'Vfi\n\ujdx 

V J ijJ J 


= / q*{v + w){<p + f)) dx - / V{DVv\^)(l)\^ dx 


V{DVw\n\^)fi\n\^ dx 


' fl\uj 

X^[V(DVn|a;)] - Xov[V(DVw|q\^)] + q *{v + w)^ {fi + ^jJ) dx 

< max{l, ||g*|U-(o)}(||t^||p + lklly)||0 + ' 011 ^ 2 ( 0 ) • 

where Xco ■ -^^(<^) Xov : \ cu) —)■ denote the operators defined by 

the respective extension by zero to all of f2. 

The operator A{u) can be estimated as follows: For all u G X we get 

\\A{u* + v)\\q^x = \\{u* + v)q\\Q^x = sup 

qeQ,q^O II^IIQ 

\\u* + v\\L2(Q)\\q\\L^n „ ||w* + t'||i,2(Q)||g||j^s(n) 

S sup -rr-S sup -[]—[]- 

q&Q,q^O ||?||H“(n) 96(3,97^0 

< CH‘>^L°°{\\u*\\L2{n) + ||r’||L2(n)), 

i.e. dl]) in Assumption is satisfied with Ca = max{l,supj>o lli(^)} which 

by continuity of the embeddings V ^ X and V + V X implies ([3]) and (|5]). 

The nullspace of G and its orthogonal complement are given by 

Af{G) = {w E L^(f2) I = 0} = {zc G L^(f2) | supp(tc) C \ a;} , 

Af{G)-^ = {n G I v\n\^ = 0} = {z; G L^(n) | supp(z;) C u} , 

and the respective projections are defined by 

Ru = Xlo[u\J\ , Pu = u - Ru = xn\Lo[u\n\Lo] ■ 

We define the operators M, N and the spaces VX,VX as follows. 

{Mv,fi)x = [{\D\Xv\^fX(j)\^dx+[v(l)dx WvEV,(j)EV^ (44) 


{Nw,fi)x= / (|T>|Vzz;|n\^)'^Vz/)|n\^dx+ impdx ^wEV^fiE 

J Q\uj J ft 

(making use of the fact that the spaces V°°, 14°° are dense in 14,14, respectively). Hence 

assuming that D does not change its sign on uj and on H \ a; (D > 0 a.e. on a; or D < 

0 a.e. on u) and (T> > 0 a.e. on H \ a; or D < 0 a.e. on H \ a;) we get 


II nj II {Mv,(l))x ^ I, I 

||Mn||x = sup -x—r. -< ItI 




II0IIl2(o) 


V 


VveV 


II Ar II {Mw,ij)x ^ I 

||Azz;||x = sup . .< llzcl 


V'GV"®, V’7^0 


II'*AIU2(o) 


V 


VweV 
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(Mn, v)x = II v^VnU||i2(^) + |Mli2(o) =: lk||~ 'iveV 

{Nw,w)x = |lv^Vu)|nv||i2(ov) + \\w\\l2^n) =■ Ikife VweV 

Since in this case G has closed range we need not regnlarize in case of noisy data, i.e., 

we can set a = 0: 




= = Xuj[z^] , 4 = - z) = - z] 



In our implementation we consider the one dimensional case with domain 12 = (0,1). 
The right hand side / is given by f{t,x) = + g*(x)) sin(7ra:), where the 

exact parameter q* is a quadratic polynomial, q* = 0.025x^ — 0.025x. For simplicity the 
diffusion coefficient is chosen to be constant, D = 1. 

For solving the partial differential equation system ([6])-([8]) we derive its variational 
formulation and discretize the spaces Q and X by cubic Hermite basis functions and 
fjj for j = 2,..., — 1 on a uniform mesh 0 = Xi < X 2 < ... < Xx = 1, where = 31 in 

our case. 


(j)j{x) 



if X G {xj-i, Xj) 


0 else 


fifix) 



if X G (xj_i, Xj) 



if X G (xj,Xj+i) 


The reason for using such high order spaces is the required regularity on arguments 
of the operators M, N according to fHT)) . flT5]) . After using these as ansatz and test 
function in the variational formulation for space discretiztion, we solve the resulting 
ODE System with an implicit Euler method with step size ht = 0.6. The interesting 
cases are those with partial observations and noisy data. 

In our experiments we just employ a simple heuristic choice of p and v. In case 
of partial observations we used the lower bound for /r. Assumption 13.41 where the 
constant Ci shows up. Therefore we solve the optimization problem min^j \\Ru — Ru*\\‘j^, 
where the constant ci varies in decimal steps between 0.1 and 1000, in order to hnd an 
appropriate /r. For the case of noisy data the constants p and z/ vary between 0.1 and 
1000 and we solve the minimization problem \\Ru — Ru*\\\. This approach will 

be enhanced in future work. 

To investigate the case of partial observations we restrict the data to the subinterval 
u = (0.3,0.87) of D = (0,1). The results for this case are shown in Figure ITTl There 
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exact and estimated solution q for different time steps t exact and estimated soiution u for different time steps t 




Figure 1. The parameter estimate q{t,x) (left) and the state estimate u{t,x) (right) 
at times t = 0,6,15,30,45,60. 


on the left the estimated parameter for different times varying from [0,60], are shown 
starting with g(0) = 0. The estimated parameters are the lines with markers, whereas 
the straight line indicates the exact parameter. On the right, the estimated (lines with 
markers) and exact (straight lines) state for different time steps are displayed. Althongh 
we have jnst partial observations, the state is estimated qnite well also in the nnobserved 
region O \ w. One can see that also the estimated parameter gets close to the exact 
one, bnt it is shifted to the right, which is dne to the fact, that data are just given 
on the nonsymmetric interval u. Note that we do not know whether the persistence of 
excitation condition is satished here, which is in fact hard to verify in general. 

For the noisy data case we assumed to have data with Gaussian noise with different 
noise levels 6 = 1%, 5%, 10%. In this case of irregular noise, according to section ??, 
parameter and state convergence cannot be proven if T* < oo, so we expect closeness 
only for times satisfying condition (133|) . This can be seen in the numerical results as 
well, because the error is increasing from a certain time instance on, which corresponds 
to the semiconvergence phenomenon in regularization. As one might expect the time 
where the error starts to grow again gets smaller as the noise level increases. In Figure 
14.11 the above row shows the estimated parameter for the three different noise levels 
6 = 1%, 5%, 10% for different times up to f = 75. 

In Figure im the lower row displays the errors of the estimated observed state and 
parameter (||e||Q + ||r||^) as in Proposition 13.131 developing over time. For small noise 
6 = 1% the estimated parameter gets close to the exact one, and also the error decreases, 
as time proceeds, whereas for larger noise 5 = 5%, 10% the estimated parameter hrst 
gets close to the exact parameter up to a certain time instance and then it drifts away 
again, hence the error increases. 
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Figure 2. Estimated parameter q{t,x) at different times and the observed and 
parameter error (||e||Q + ||r||3(-) of Proposition 13.131 for different noise levels S = 
1,5,10%. 

5. Conclusion and Outlook 

In this paper we have developed and analyzed an online parameter identification method 
for time dependent problems. The main idea was to formulate an alternative dynamic 
update law for the state and an additional one for the parameter estimate. We showed 
that the solution of this alternative system of differential equations is well defined 
and that it converges to the exact parameter and state. The proofs were done for 
the case of exact data as well as for the case of noisy data and smooth noisy data. 
The main advantages of this method are, that it imposes less restrictions on the 
underlying model compared to existing methods and that it is also applicable in case of 
partial observations. In a numerical example we showed the performance of our online 
parameter identification method. Here the results could be improved by finding optimal 
values for /i and u. 

Another future goal is to consider time dependent parameters. This would mean 
that the model itself contains a dynamical update law for the parameter and therefore 
the online parameter identification method has to be adapted, so that this is taken into 
account. 
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